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Abstract 

The energetic beam of (spin and isospin zero) a-particles remains a very ef- 
ficient probe for the nuclear isoscalar giant resonances. In the present work, 
a microscopic folding model study of the isoscalar giant resonances in 208 Pb 
>• ' induced by inelastic a+ 208 Pb scattering at = 240 and 386 MeV has been 

oo : 

■ performed using the (complex) CDM3Y6 interaction and nuclear transition 
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densities given by both the collective model and Random Phase Approxi- 
mation (RPA) approach. The fractions of energy weighted sum rule around 
the main peaks of the isoscalar monopole, dipole and quadrupole giant res- 
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onances were probed in the Distorted Wave Born Approximation analysis 



of inelastic a+ 208 Pb scattering using the double-folded form factors given 
by different choices of the nuclear transition densities. The energy distribu- 
tion of the EO, El and E2 strengths given by the multipole decomposition 
analyses of the (a, a') data under study are compared with those predicted 
by the RPA calculation. 
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1. Introduction 

Isoscalar giant resonances jl| in medium and heavy nuclei are the pro- 
nounced manifestation of nuclear collective motion and, hence, they carry 
important information about the dynamics of the nuclear excitation process 
and the properties of the nuclear Hamiltonian. Although their systematic 
study started more than three decades ago, a number of challenging ques- 
tions still remain open. For example, the observed compressional modes, like 
the L = isoscalar giant monopole resonance (ISGMR) or the L = 1 isoscalar 
giant dipole resonance (ISGDR), provide the optimal route to determine the 
nuclear matter incompressibility K^, a key quantity specifying the equation 
of state of nuclear matter. However, the compatibility of the values de- 
duced from these two collective compressional modes is still under discussion 
(see, e.g., Refs. and references therein). The isoscalar giant quadrupole 
resonance (ISGQR) and the isovector giant dipole resonance (IVGDR) are 
well known in stable nuclei but it is remains unclear how these two modes 
evolve in the unstable neutron-rich nuclei. 

Since the first observation of ISGMR in 208 Pb in the 70's of the last cen- 
tury 0, [5, (|, this compressional excitation mode of 208 Pb has been investi- 
gated in numerous experimental and theoretical studies due to its fundamen- 
tal importance in the determination of the nuclear matter incompressibility 
Kqo. Among the experimental tools, the spin and isospin zero a-particle re- 
mains the best probe for the isoscalar giant resonances with AS = AT = 
and the most pronounced observations of the ISGMR in 208 Pb have been 
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made so far in the (a, a') experiments, like the recent high-precision mea- 
surements of inelastic a+ 208 Pb scattering at E^b = 240 MeV by the Texas 
A&M University group [7| and 386 MeV by the Osaka group 0,3- The 
exploratory theoretical studies of the isoscalar dipole mode were done in the 
early 80s by M.N. Harakeh and A.E.L. Dieperink [lOj, as well as by N. Van 
Giai and H. Sagawa 11], but the first direct observation of the ISGDR in 
the 208 Pb(a, a') reaction has been made only years later by Davis et al. 12]. 
Like the ISGMR, the knowledge about the ISGDR is of vital importance for 
the determination of the nuclear incompressibility ifooj^, 3]. Therefore, both 
the (a, a') experiments at 240 MeV \f\ and 386 MeV 8, J were aimed at an 
accurate measurement of the ISGDR strength distribution in 208 Pb. One of 
the main problems in the experimental study of giant resonances has been, 
and still is, the difficulty to disentangle different modes when their energies 
overlap. For example, the excitation energy of ISGMR in 208 Pb has been 
accurately determined from the most forward part of the (a, a') cross section 
measured at 240 MeV to be E « 13.96 ± 0.20 MeV Q while the (a, a') ex- 
periment at 386 MeV has deduced E m 13.4 ± 0.2 MeV {^J using essentially 
the same method of multipole decomposition analysis (MDA) to disentangle 
the ISGMR peak from a mixed spectrum of different (AS = AT = 0) excita- 
tions including, in particular, the low-energy peak of ISGDR at an excitation 
energy around 13 MeV 

Moreover, there exists a basic problem when experimental results need 
to be compared with those predicted by a theoretical structure model. It is 
common, in the theoretical works, to calculate the strength function S(E) 
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associated with a given nuclear transition operator Q, that is, 

S(E) = ^26(E-E U ) | <u\Q\0> | 2 , (1) 

V 

where v labels a complete set of final states which can be excited by acting 
Q upon the ground state |0 >. In terms of single-particle degrees of freedom, 
the isoscalar L-multipole operator is given by 

A 

QLM = Y. r ^{h). (2) 

1=1 

In the case of the isoscalar monopole and dipole modes one has to replace rf 
by rj" +2 . The predicted strength function ([1]) is then used to compare with 
the corresponding experimental strength distribution deduced for the con- 
sidered L-multipole isoscalar excitation. In reality, however, the measured 
inelastic scattering data are inclusive spectra over a wide energy range which 
contain the strengths of isoscalar excitations with different multipolarities as 
well as the contamination by the continuum background. The 'experimental' 
strength function for a given L-multipole excitation depends strongly, there- 
fore, on both the method to exclude continuum background and the MDA 
to disentangle contribution of a given multipole from the inelastic scattering 
spectrum. As a result, the comparison between the theoretical and experi- 
mental strength functions can be made in many cases only qualitatively. 

Nevertheless, a more direct and quantitative comparison between theory 
and experiment is possible if the nuclear structure information is accurately 
included into a microscopic description of the inelastic (a, a') scattering cross 
section measured for a given peak of the resonance energy, within either the 
distorted-wave Born approximation (DWBA) or coupled-channel formalism 
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[ljj. The well-known portal of this procedure is the double-folding model 
(DFM) which uses the nuclear ground state and transition densities of the 
a-projectile and target nucleus and an effective nucleon-nucleon (NN) inter- 
action to calculate the a-nucleus optical potential (OP) and inelastic scatter- 



ing form factors 



T) of different multipolarities for the DWBA calculation 



see, e.g., Ref. 141 ] and references therein). We note that the single- folding 



method 



151 ] has been widely used to calculate the a-nucleus inelastic scat- 



tering FF, using an appropriate a-nucleon interaction and nuclear transition 
densities given by the collective model, for the multipole decomposition anal- 
ysis of the experimental spectrum within DWBA {?], s[ 0]. In this respect, 
the DFM calculation of the a-nucleus inelastic scattering FF using a real- 
istic effective NN interaction can be used to probe the transition strength 
extracted from the MDA for a given isoscalar excitation as well as the re- 
liability of different choices for the nuclear transition densities. The latter 
aspect is quite essential because the collective model transition density has 
been shown to be reasonable for the collective modes which rule, 
concentrated in a limited energy region while other parts of the spectrum 
could be dominated by pure particle-hole states or states with intermediate 
character. 

We note further that the existing structure models have been substan- 
tially improved in recent years and the nuclear linear response theory has 
been used successfully to describe the excitation of vibrational modes. In a 
magic nucleus like 208 Pb where the pairing does not manifest itself, this linear 
response approach is also known as Random Phase Approximation (RPA) 
method. RPA has been formulated many years ago, but only recently the 
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fully selfconsistent calculations without crude approximations have become 
available. Therefore, it is also timely to study the inelastic a-nucleus scat- 
tering data measured recently for the known isoscalar modes to see how the 
RPA nuclear wave functions can be probed in the folding + DWBA analysis 
of these data. To this goal, we have performed in the present work a detailed 
folding model analysis of the high-precision inelastic a+ 208 Pb scattering data 
measured at = 240 MeV 7] and 386 MeV |g[ 9L using a complex version 
of the density dependent CDM3Y6 interaction 16j and nuclear transition 
densities given by both the collective model (CM) and the RPA calculation. 
After a brief overview of the theoretical formalism in Sec. [21 results of the 
DFM + DWBA analysis of the considered (a, a') data using the CM transi- 
tion densities are presented and discussed in Sec. El The RPA description of 
the E0, El and E2 strength distributions and DFM + DWBA results given 
by the RPA transition densities are discussed in Sec. HI The main conclusions 
and perspectives are given in the Summary. 

2. Formalism 

In this Section, we briefly describe the theoretical model used to calcu- 
late the inelastic (a, a') cross sections in the DWBA. As mentioned above, 
our microscopic study of the inelastic a+ 208 Pb scattering is based on the 
double-folding model [3] which uses the nuclear ground state and transition 
densities of the a-projectile and target nucleus and an appropriate effective 
NN interaction to calculate the a-nucleus OP and inelastic scattering FF 
for the DWBA analysis. The nuclear structure information on the isoscalar 
giant resonances in 208 Pb is embedded in the nuclear transition densities used 
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in the folding calculation of the inelastic scattering FF. Two choices of the 
nuclear transition densities were used in the present work: the phenomeno- 
logical transition densities given by the collective model (discussed below in 
Sec. 12.3ft and microscopic transition densities given by the RPA approach 
(discussed in Sec. 



2.1. Effective density dependent N N interaction 

Among various choices of the effective NN interaction, a density depen- 
dent version of the M3Y interaction (dubbed as CDM3Y6 interaction [161 ]) 
has been used successfully in the folding model analyses of the (refractive) 
elastic and inelastic a- nucleus scattering (see the recent review in Ref. {itI]). 
The density dependent parameters of the CDM3Y6 interaction were carefully 
adjusted in the Hartree-Fock (HF) scheme to reproduce the saturation prop- 
erties of nuclear matter [l6]. The first version of the CDM3Y6 interaction 
is real and can be used to predict the real OP and inelastic scattering FF 
only. To avoid a phenomenological choice of the imaginary OP and inelastic 
scattering FF, we have supplemented the real CDM3Y6 interaction with a 
realistic imaginary density dependence whose parameters were determined 
based on the Brueckner Hartree-Fock (BHF) results for the nucleon OP in 

n 

nuclear matter by Jeukenne, Lejeune and Mahaux (JLM) [18(. It has been 



shown in our recent work |l9j that the same form of the CDM3Y functional 
\\\ can be used to obtain the density dependence of the imaginary term. 
Thus, the complex CDM3Y6 interaction used in the present folding model 
analysis is determined as 

Re(Im) w D (ex) (E, p, s) = F R(I) (E, p)v d(ex)(s) , (3) 
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F X (E, p) =C X [1 + a x exp(-/3 x p) - j x p], x = R, I. (4) 

The radial parts of the direct and exchange interactions Ud(ex)(s) were kept 

14J from the M3Y interaction based on the G-matrix 



unchanged, as derived 



elements of the Paris NN interaction [20], in terms of three Yukawas 

Ms) = 11061.625 eXp( " 4s) - 2537.5 eXp( - 2 - 5g) , (5) 

4s 2.5s 

exp(-4s) exp(-2.5s) exp(-0. 7072s) 

Vex ( s ) = -1524.25—^ ^-518.75—^ - - 7.8474 1 K - - '-. 

v ; 4s 2.5s 0.7072s 

While parameters of the real density de pen dence Fr were taken from the 



original HF calculation of nuclear matter 



161 ] , those of the imaginary density 



dependence F\ were adjusted iteratively until the HF result for the imaginary 



nucleon OP in nuclear matter agrees reasonably with the JLM result [18| as 
well as the shape of imaginary folded OP becomes close to the phenomeno- 
logical Woods-Saxon imaginary OP found at each energy. All parameters 
of the complex density dependence are given in Table 1. We note that the 
dynamic change in the density dependence F x (p) caused by the excitation 
of the target is taken into account properly in the folding calculation using 
method given in Ref. [l^ . 

2.2. Double- folding model 

The generalized double-folding model of Ref. [14| was used to evaluate 
the complex a-nucleus OP and inelastic scattering FF from the following 
HF-type matrix elements of the CDM3Y6 interaction (jSJ)-© between the 
projectile nucleon i and target nucleon j 

U a ^a*= ^2 [< u'I^dIu > + < u'I^exI^* >], (6) 

i£a;j£AJ'£A* 
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where A and A* are states of the target in the entrance- and exit channel 
of the a-nucleus scattering, respectively. Thus, Eq. (jSJ) gives the (diagonal) 
elastic OP if A* = A and (nondiagonal) inelastic scattering FF if otherwise. 
The (local) direct term is readily evaluated by the standard double-folding 
integration 

U D (E,R) = J p a (r a )p A (r A )v D (E,p,s)d 3 r a d 3 r A , 

s = r A -r a + R. (7) 

The antisymmetrization gives rise to the exchange term in Eq. ([6]) which is, 
in g eneral, nonlocal in the coordinate space. However, it has been shown 
flX \v\ that an accurate local equivalent exchange potential can be obtained 



using the local WKB approximation 13| for the change in relative motion 



induced by the exchange of spatial coordinates of each interacting nucleon 
pair 



U E x(E,R) = J p a (r a ,r a + s)p A (r A ,r A - s)v E x(E,p,s) 

fiK(R)s\ ,o ,o 

x exp I y M ' J d 3 r Q d 3 r A . (8) 

Here K(R) is the local momentum of relative motion determined from 

K\R) = j^[E c . m . - Re U (E,R) - V C (R)\, (9) 

where p is the reduced mass, M = 4v4/(4 + A), E c _ m _ is the scattering energy 
in the center-of-mass (cm.) frame, Uq(E,R) and Vc(R) are the nuclear 
and Coulomb parts of the real a-nucleus OP, respectively. The calculation 
of Untyx ) is done iteratively based on a density-matrix expansion method 
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23j |. All technical details of the folding calculation of C/d(ex) are the same 



as those given in Ref. [14J], excepting the use of a realistic local approximation 
for the transition density matrix suggested by Love [24| and a recoil correction 



to the exchange term (jSJ) suggested by Carstoiu and Lassaut 



25]. 



2.3. Collective model for the nuclear transition densities 

To calculate consistently both the OP and inelastic scattering FF for the 
a+ 208 Pb system one needs to represent the target density in terms of the 
ground state (g.s.) and transition parts as p(r) = Po(r) +5p(r). The explicit 
expression of the inelastic scattering FF for a given isoscalar excitation (see 



Ref. 



141 ]) can be deduced from the double-folding integrals fl7D-flSJ), using 



the following multipole decomposition of Sp(r) 

5p{r) = Y,C L 5p L (r)\i L Y LM (r)Y 



(10) 



LM 



where Cq = yAn and Cl—1 for L 7^ 0. Given the strong collective nature of 
the isoscalar giant resonances, macroscopic methods to construct the nuclear 
transition density of the 2 L -pole isoscalar excitation 5pi(r), based on the 
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26] 



collective model, are widely used in the folding model calculation 

and multipole decomposition analysis q y, 0, h 3 of the (a, a') data. 

For the isoscalar giant resonances with L > 2, we adopt the so-called 
Bohr-Mottelson prescription [^tJ to construct the transition densities 

dp (r) 



SpL(r) = -S L - 



dr 
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Here po( r ) is the g.s. density and 5l is the deformation length of the con- 



sidered isoscalar excitation, 
distribution with parameters 



he g.s. density of 208 Pb was taken ClS db Fermi 



281 ] chosen to reproduce the shell-model density 



n 



for 208 Pb. Within the isoscalar assumption 26] the same deformation 
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length 5l is employed, as a rule, for both the neutron and proton parts of 
the nuclear transition density pip . For the low- lying excitations, like the 
first 3~ state in 208 Pb considered below, the deformation length is normally 
determined [14J from the measured electric transition strength B(EL). In 
terms of the energy weighted sum rule (EWSR) for the operator (J2J), if a 
single state \u > at the excitation energy E v exhausts 100% of the isoscalar 



EWSR then the corresponding deformation length is determined 



as 



h 2 4tt L(2L + If < r 2L ~ 2 > 
2mAE u (L + 2) 2 < r L ~ l > 2 



< r 



L-l 



> 



J>o< 



r r 



L+l 



dr 



J po{r)r 2 dr 



:i2) 



In the case of ISGMR, the pure breathing mode (or scaling) assumption 29] 
is used to construct the nuclear transition density 

dp (r 



-<$n 



3po(^j 



dr 



(13) 



If an isoscalar monopole state \u > at the excitation energy E u exhausts 
100% of the monopole EWSR then its deformation length is determined 26] 

as 

n9/ . h 2 An 1 . , 

5 2 (E U ) = — — — ^ . (14) 



2m AE V < r 2 >' 

Another special case is that of the isoscalar dipole excitation for which 
a macroscopic model based on the compressional hypothesis, with a proper 
center-of-mass subtraction, has been suggested by Harakeh and Dieperink 
flo| . Dropping the high-order term e which is negligible for A > 20 10], the 
transition density of an isoscalar dipole state is written as 



3r 2 4- + I0r--<r 2 >4- 
dr 3 dr 



Po{r) 



(15) 



where R is the half-density radius of the g.s. density distribution po( r )- If 
an isoscalar dipole state \v > at the excitation energy E v exhausts 100% of 
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the 



id, 



WSR for the dipole operator, with spurious cm. oscillation subtracted 



111 ], then its deformation length is 



n ; mAE v 



it 4 2 2 

11 < r 4 > < r > z 



(16) 



The CM transition densities (ITT|) - (IT6"|) are normalized in our calculation to 
describe the excitation process |g.s. >^ \v >, that is, they correspond to 
the upward transition amplitude. In this way, the corresponding isoscalar 
transition strength is Sl = \Ml\ 2 with the transition moment determined as 



M L = J dr r L+2 5p L (r) if L > 2, 
M L = J dr r A 5p L {r) if L = 0, 
M L = J dr(r 3 -^<r 2 > r^j r 2 5p L (r) if L = 1. (17) 

2-4- Microscopic RPA transition densities 

Despite a certain success of the collective model transition densities in 
numerous folding model studies of the isoscalar giant resonances induced by 
inelastic a-nucleus scattering, there is no firm experimental evidence vali- 
dating their use. In the case of ISGMR, for example, there are only some 
results of structure calculations showing that ([TBI is a good representation 



of the ISGMR transition density in the surface region |30(. Moreover, the 
radial shapes of the collective model transition densities are assumed to be 
independent of the excitation energies, which is surely not the case in the re- 
ality. Therefore, the folding model analysis based on the transition densities 
(for a given isoscalar mode) calculated selfconsistently at different excitation 
energies by a microscopic RPA or quasiparticle RPA (QRPA) approach is 
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expected to provide a complementary and useful insight. While the QRPA 
transition densities have been shown to give reasonably good results in many 
cases like, e.g., in the folding model study of the lowest 2 + states in the Sul- 



fur isotopes induced by inelastic proton scattering [3l|, QRPA or RPA do 
not systematically provide good results for the low-lying isoscalar excitations 
when strong anharmonic effects are present. In this sense, RPA is expected 
to be more suitable for giant resonances and it is, therefore, of interest to 
probe the RPA transition densities in the present study. We note in this 
context a similar attempt done recently to study the charge exchange ( 3 He,t) 
reaction |32] . In general, the full coupling between the microscopic structure 
and reaction models should be, in our opinion, pursued more extensively. 
In the present calculations, we have chosen the parametrization set SLy5 



331 ] of the Skyrme interaction for the RPA calculation of the isoscalar states 



in 208 Pb. We first solve the HF equations in the coordinate space to con- 
struct the single-particle basis. All the radial integrals are computed up to 
a maximum radius of 22.5 fm, using a mesh of 0.15 fm. The unoccupied 
single-particle states, including those at positive energies, are obtained by 
putting the system in a large box of 22.5 fm, that is, the continuum is dis- 
cretized. A basis of particle-hole (ph) configurations is then built upon all 
occupied states, as well as the lowest unoccupied states with increasing val- 
ues of the principal quantum number n, for each allowed value of (/, j). The 
RPA matrix equations are then solved in this basis, which has been checked 
to be large enough to ensure that the appropriate sum rules are satisfied. 



The procedure has already been explained in Ref. 



34|- 



From the solutions of the RPA equations, the energies E v of the excited 
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states \u > as well as their wave functions are readily obtained. The radial 
transition density 5pL„(r) associated with a given (2 L -pole) RPA state \v > 
is given by 

8p$(r) = {K^ + Y ih U) ) < P\\YL\\h > R P (r)R h (r), (18) 

phGq 

where X and Y are the forward and backward RPA amplitudes and R(r) 
labels the radial part of the single-particle wave function. The proton and 
neutron parts (labelled by q = p, n) of the transition density flTSJ) are com- 
puted separately. The isoscalar transition strength of the RPA state \v > is 
evaluated as Sl v = \Ml u \ 2 , where the transition moment Ml u is 



M. 



Lv 



M, 



Lv 



dr r L+2 
dr r 



M, 



Lv 



dr ( r 3 < r 2 > r ) r 2 

3 



"*pg(r)+*pS(r) 
5p { L°!(r)+5p^(r) 



if L > 2, 
if L = 0, 
if L— 1(19) 



3. Double-folding model + DWBA analysis using the collective 
model transition densities 

To generate realistic distorted waves for the DFM + DWBA study of the 
isoscalar giant resonances, we first used the nuclear g.s. densities of 4 He and 

respectively, to calculate the complex 



208 



Pb taken from Refs. Ij33l and 



folded OP for the optical model (OM) analysis of the elastic a+ 208 Pb scat- 
tering data at E Uh = 240 MeV [z] and 386 MeV fl. To fine tune the complex 
strength of the CDM3Y6 interaction ([3]), renormalization coefficients Nr and 
Nj of the real and imaginary elastic folded potentials (jSJ) were adjusted by 
the OM fit to the elastic data at each energy (see OM results shown in upper 
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panel of Fig. [T]). One can see from Tabled] that the best-fit coefficient is 
rather close to unity. The best-fit N\ of about 1.4 is reasonable because the 
imaginary strength of the CDM3Y6 interaction was tuned to the BHF results 
for nuclear matter and gives, therefore, only the "volume" absorption. To ef- 
fectively account for the surface absorption caused by inelastic scattering and 
transfer reactions, an enhanced Ni coefficient is naturally needed (compare 
the dash and solid curves in upper part of Fig. [p. Our OM calculation also 
predicted the total reaction cross sections ctr very close to the experimental 
values measured at the nearby energies. Thus, the elastic distorted waves 
given by the present DFM calculation should be accurate for the DWBA 
analysis of inelastic a+ 208 Pb scattering. 

For the inelastic scattering form factor, a standard method used so far in 
the DFM + DWBA aI1 a ly ses of initio .-nude, scatter fl Q is to 
scale the real and imaginary inelastic folded FF by the same renormalization 
coefficients and Ni as those deduced from OM analysis of elastic scattering 
data. We show in lower panel of Fig. [TJ as illustration, the DWBA description 
of inelastic a+ 208 Pb scattering data for 3j~ state of 208 Pb given by the inelastic 
folded FF scaled by the same coefficients and Nj as those given in Table [TJ 
By using a deformation length Sl of the CM transition density ffTTl) chosen 
to reproduce the measured transition rate B exp (E3) ~ 611 x 10 3 e 2 fm 6 



a very satisfactory description of the inelastic scattering data for 3^ state 
of 208 Pb has been obtained. About the same good DWBA description was 
also obtained with the microscopic nuclear transition density f|T8|) given by 
the RPA calculation, without any ad hoc adjustment. We note that the 
Coulomb part of inelastic scattering FF is obtained in the present work by 
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double folding the proton parts of the 4 He g.s. density and 208 Pb transition 
density with the Coulomb interaction, using a folding method similar to that 
used for the nuclear part. 

3.1. MDA and deformation lengths for the CM transition densities 

Before discussing the DFM + DWBA results for the isoscalar giant res- 
onances we briefly recall here how the experimental transition strengths are 
determined from the multipole decomposition analysis of the measured (a, a') 
angular distributions [7, 0, At a given energy bin, the measured (dou- 
ble) differential cross section is expressed within the MDA as a superposition 
of the angular distributions calculated for different transferred angular mo- 
menta L as 



d 2 a 



(Oc.m. ; E x 



cxp. 



d 2 a 



dVLdE 



(Oc.m.; E x 



calc. 



(20) 



L=0 

Here [d 2 a / dVtdE\ c ^ Xc ' is calculated within the DWBA using the inelastic scat 
tering FF generated from the appropriate CM nuclear transition density (TTT|) 



f|T6|) by a single-folding method 15j. The CM nuclear transition densities 
entering the MDA are first determined with 100% exhaustion of the corre- 
sponding EWSR (see Sec. I2.3p . then a least-x 2 -fit procedure determines all 
ai{E x ) coefficients for the considered experimental energy bin. As a result, 
each best-fit ai{E x ) coefficient represents the fraction of EWSR exhausted by 
the corresponding isoscalar 2 L -pole excitation mode in the energy bin under 
the MDA analysis. In terms of deformation length 5l for a given excitation 
mode in the considered energy bin ol{E x ) = (5_L/5™ ax ) 2 , where <5™ ax is the 
maximum deformation length determined [see Eqs. (j!2p . f[T4l) and f|T6|) ] to 
exhaust 100% of the corresponding EWSR. 
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The MDA analysis of the inelastic a+ Pb scattering data at E^b = 240 
MeV measured by the Texas A&M University group 7| was done in the 
energy bins of 640 or 800 keV width to deduce the isoscalar EL strength 
distributions over a wide range of excitation energy. The MDA of the 240 
MeV data shows, in particular, that the full exhaustion (around 100%) of the 
isoscalar EWSR has been observed for the ISGMR and ISGQR. The main 
ISGMR peak has been accurately determined from the 240 MeV data to be at 
E x « 13.96±0.20 MeV with a width T « 2.88±0.20 MeV, and fragmentation 
of the E0 strength up to about 20 MeV has been observed. The MDA analysis 
of high-precision inelastic a+ 208 Pb scattering data measured by the Osaka 
group at £i a t, = 386 MeV js, 9] (done in the energy bins of 1 MeV width) 
has shown a much stronger fragmentation of the E0 strength over excitation 
energies well above 30 MeV, and a less pronounced ISGMR peak (observed 
at E x « 13.4 ± 0.2 MeV with a wider width r « 4.0 ± 0.4 MeV). 

One could reproduce the ISGMR peaks observed in the (a, a') experiment 



at 240 MeV |7f and 386 MeV [8j, |9j in the microscopic structure models, using 
either nonrelativistic or relativistic functionals which give Koo « 240 ± 20 
381 ] . This is why some consensus has been reached on this 



MeV 



33, 



empirical value for K^, where the error of ± 20 MeV is not simply associated 
with the experimental uncertainty on the ISGMR energy, but rather with 
our still incomplete understanding of the structure of energy functionals (in 
particular, of their density dependence). In this connection, we note that a 
pure experimental discrepancy of 500 keV in the observed ISGMR peaks could 
also result in a difference of AK^ « 20 MeV in the values extracted from 
empirical formulas relating the ISGMR peak in 208 Pb and nuclear matter 
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incompressibility K 



0, y 



41] . Such a difference is quite sizeable and 
hinders any further theoretical modelling of the energy functionals. 

If one uses a microscopic structure approach to determine from the 
ISGMR data, the location of -Eisgmr will affect the deduced value. In 
about the same waythe observed ISGDR peak may also be directly related 
to the value 1, 2, jj. Given such a vital importance of the ISGMR and 
ISGDR excitations in determining the nuclear matter incompressibility 
and the fact that a simple single-folding method [15J was used to calculate 
the ct+ 208 Pb inelastic scattering FF for the MDA analyses of Refs. % |8|, |9|, 
we deem it necessary to probe the ISGMR and ISGDR strength distributions 
extracted from these two experiments again in our DFM + DWBA approach. 

It is complementary to note that the density dependent CDM3Y6 in- 
teraction used in the present DFM calculation was parametrized in the HF 
scheme to reproduce = 252 MeV at the saturation density of symmet- 



ric nuclear matter 
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171 ]. and it has been successfully used in numerous 
OM and DWBA analyses of elastic and inelastic a-nucleus scattering. The 
isoscalar EL strengths (in terms of exhausted fractions of the correspond- 
ing EWSR for L = 0, 1, 2, 3) given by the MDA analyses of Refs. Q, Q 
for the main peaks of the ISGMR, ISGQR and ISGDR together with those 
predicted by the RPA calculation are presented in Table [2j Each energy bin 
in Table [2] has been chosen so that the strongest EL strengths around the 
main resonance peaks deduced from the two experiments can be used con- 



sistently in the same DFM + DWBA calculation. For example, the 
peak was found by the MDA of 240 MeV data [7J and 386 MeV data 



SGMR 



9] 



m 



the energy bins centered at E x 14.1 and 13.5 MeV, respectively, and the 
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corresponding EO strengths should be studied in the same DFM + DWBA 
analysis. The full energy distributions of the EL strengths for L — 0, 1, 2 are 
shown below in Sect . HI 



3.2. Isoscalar EL strengths near the ISGMR peak 

In Fig.[2]the inelastic a+ 208 Pb scattering cross section at E^ = 240 MeV 
measured for the 640 keV-wide energy bin centered at E x = 14.1 MeV [3] 
are compared with the DFM + DWBA results given by the collective model 
transition density (1T31 - (1T4"I) . In this energy bin, the EO strength deduced 
from the 240 MeV data is strongest and exhausts about 37.6% of the EO 
EWSR. The isoscalar El strength is quite significant (6.3% of the El EWSR) 
in this energy bin due to the dipole strength coming from the low-energy 
peak of ISGDR located around 13 MeV and affecting significantly the total 
angular distribution. The isoscalar E2 strength of about 6.6% of the E2 
EWSR as well as no contribution from isoscalar E3 excitation were found 
in this energy bin [7j|. One can see in the lower panel of Fig. [2] that the 
measured inelastic scattering cross section is reasonably described by the 
DFM + DWBA calculation using the CM transition densities scaled to the 



0. 



Given the 



isoscalar (L = 0, 1, 2) strengths given by the MDA of Ref. 
monopole and dipole angular distributions oscillating out of phase, a smooth 
angular distribution seen in the 240 MeV data shows clearly the mixture 
of the El strength from the low-energy peak of ISGDR in the considered 
energy bin. Since the MDA of the Osaka data js, ^ has given the strongest 
E0 strength in the 1 MeV energy bin centered at E x = 13.5 MeV, we found 
it appropriate to use the isoscalar EL strengths deduced for this energy bin 
to construct the CM transition densities for our DFM + DWBA analysis of 
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9j much stronger (up to 



the 240 MeV data for E x = 14.1 MeV. According to the MDA of the Osaka 
data, the E0 strength is strongly fragmented over a wide energy range and 
only about 16% of the E0 EWSR has been located in the energy bin around 
the ISGMR peak at E x = 13.5 MeV. While the isoscalar dipole strength 
deduced from the 386 MeV data is quite close tojthat deduced from the 240 
MeV data, the isoscalar E2 strength was found 
15% of the E2 EWSR) in the energy bin around E x = 13.5 MeV. On top 
of that, about 3% of the isoscalar E3 EWSR strength was also observed by 
Uchida et al. j], ^] in this energy bin. The DFM + DWBA description of the 
240 MeV data for E x = 14.1 MeV given by the CM transition densities scaled 
to the isoscalar (0 ^ L ^ 3) strengths taken from Refs. {§], 9] are shown in 
upper panel of Fig. [2J Although a weaker E0 strength deduced by Uchida 
et al. 8j, |9J gives a monopole cross section more than 2 times smaller than 
that obtained with the E0 strength deduced by Youngblood et al. [7|], the 
overall description of the 240 MeV data given by the isoscalar EL transition 

n q 

strengths taken from Refs. [8J, |9[ remains satisfactory due to stronger E2 and 
E3 contributions. However, the lack of the £"0 strength can still be seen in 
the DWBA description of data points at the most forward angles shown in 
the upper panel of Fig. [2J 



It is natural also to check 
aka data at E^ = 386 MeV 



the DFM + DWBA description of the Os- 



39) based on the same EL transition 



strengths as discussed above. In Fig. [3] the inelastic a+ 208 Pb scattering 
cross section at £j a b = 386 MeV measured for the 1 MeV-wide energy bin 
centered at 13.5 MeV are compared with the DFM + DWBA results given 
by the CM transition densities (fl3"j) - (JTlj) built upon the same isoscalar EL 
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transition strengths as those used in Fig. [2J Except for the two data points 
at the most forward angles which are fairly described by the EL strengths 
given by the MDA of the 240 MeV data [7], the present DFM + DWBA 
results strongly underestimate the measured data over the whole angular 
range. About the same picture has been found for the energy bin centered 
at E x = 14.5 MeV when the DFM + DWBA results obtained with the EL 
strengths deduced for this bin are compared with the inelastic a+ 208 Pb scat- 
tering data at E^b = 386 MeV |39|. Such a big gap between the calculated 
and measured cross sections seen in Fig. [3] is not unexpected because the 
contributions by the excitation modes of higher multipoles (L > 3) are not 
taken into account in our DFM + DWBA calculation. We recall that the 
authors of Refs. [s, 39] were able to measured the (a, a') energy spectrum 
for the lead target without any contamination from the instrumental back- 
ground by using the high-resolution magnetic spectrometer Grand Raiden, 
and the MDA analysis of the 386 MeV data has been done for all multipoles 



up to L max = 14 



39] . Since the number of fitting parameters is quite 



large in this case, it is not excluded that some continuum background com- 
ing from other quasi-elastic processes like the pickup/breakup reactions has 
been simply approximated by the high-multipole terms in the MDA series. 
Moreover, any strong particle-hole EL excitation with 9 ^ L ^ 14 in the 
energy region around 13-14 MeV would be unlikely from the structure point 
of view. In contrast to the Osaka experiment, the MDA of the 240 MeV data 
by the Texas A&M University group was done only after a broad contin- 
uum background (presumably caused by the high-multipole EL excitations 
and pickup/breakup reactions) has been substracted [7J. As a result, the 
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MDA of the 240 MeV data has been performed with less fitting parameters. 
The fact that present DFM + DWBA results describe the 240 MeV data 
reasonably using the EL strengths deduced by the Texas A&M University 
group j3] seems to indicate that the single-folding method [ij]] used in their 
MDA of inelastic a-nucleus scattering is quite reliable. However, a closer 
inspection of the lower part of Fig. [2] and Fig. 5 in Ref. \l\ shows that the 
present DFM + DWBA results are slightly underestimating the data points 
compared to the MDA results obtained with the single-folding method. To 
explore such a difference in the case of 386 MeV data, we have compared 
in Fig. H] the present DFM + DWBA results with those given by the MDA 
of the (a, a') data measured for the energy bin centered at E x = 14.5 MeV 



l39j and found that the DFM + DWBA cross sections are indeed lower 
than those given by the MDA (with the relative difference in the calculated 
total cross sections ranging from ~ 57% at the most forward angles to ~ 44% 
at 6 c . m . « 14°). In particular, the difference between the DFM + DWBA 
cross section for L = 2 and that given by the MDA is very alarming, as 
it is ~ 185% at the most forward angles and ~ 69% at 6 c . m . ~ 14° (see 
dash-dotted curves in Fig. H]) . Such a difference is quite significant and can 
result in sizable differences in the EL transition strengths deduced from the 
MDA of inelastic a-nucleus scattering using either single- or double-folding 
method. Since the present DFM approach is much more advanced compared 
to the single-folding method 15] used so far in the MDA, it is not excluded 
that the EL strengths near the ISGMR or ISGQR peaks deduced from the 
MDA are somewhat underestimated. In any case, the use of the present 
DFM approach in the MDA of future a-nucleus scattering data is strongly 
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recommended. We note further that while the ISGMR peaks observed in 
the two (a, a') measurements seem to be in a reasonable agreement with 
the existing database Qfl Q , the difference in the observed distributions of 
the isoscalar monopole strength is striking (about a factor of 2 in the most 
important energy interval). It remains, therefore, an interesting challenge to 
future experiments to confirm whether the £"0 strength is mainly localized 
at the excitation energies below 20 MeV [3] or widely fragmented to energies 



beyond 30 MeV 



3. 3. Isoscalar EL strengths near the ISGQR peak 

The ISGQR at E x » 10 ~ 11 MeV in 208 Pb is perhaps one of the most 
studied isoscalar giant resonances in nuclei. Nevertheless, like in the case of 
the ISGMR, the isoscalar E2 strength distributions observed in the (a, a') 

□ nn 

experiments at 240 MeV [7[ and 386 MeV [8|, [39J are sizeably different. For 
example, the E2 strength has been shown by the MDA of the 240 MeV data 
to be concentrated mainly near the ISGQR peak and slightly spread over 
the energies below 21 MeV, exhausting 100 ± 13% of the E2 EWSR 7|. In 
contrast, the E2 strength given by the MDA of the 386 MeV data |sj, |39| 
is broadly spread from about the same ISGQR peak to energies beyond 30 
MeV and, hence, exhausts more than 200% of the E2 EWSR (based on a 
direct integration of the tabulated E2 strength made available to us by the 
authors) . 

In Fig. 0the inelastic a+ 208 Pb scattering data at E\^ = 240 MeV mea- 
sured for the energy bin centered at E x = 10.3 MeV j3] are compared 



with t 
Refs 



re DFM + DWBA predictions based on the EL strengths taken from 



[39( | . In this energy bin, the E2 strength deduced from the 240 
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MeV data is stronger than that deduced from the 386 MeV data and exhausts 
about 20% of the E2 EWSR, and the DWBA description of the measured 
angular distribution is better (see lower panel of Fig. [5]) if the CM transition 
densities are scaled to the EL strength deduced from the 240 MeV data Q]. 
We note that a weak isovector El strength was also included into the MDA 
of the 240 MeV data to achieve a better DWBA fit to the data measured for 
the two energy bins centered at 10.3 and 14.1 MeV [7]. Given a dominant 
contribution by the isoscalar EL strengths (0 ^ L ^ 3) to the (a, a') cross 
sections (see lower panels of Figs. [2] and [5]), we have chosen not to include 
the isovector El mixing into the present DFM + DWBA calculation in or- 
der to show explicitly the role of the isoscalar excitation modes. We have 
also observed (lower panel of Fig. \5§ that the present DFM + DWBA re- 
sults slightly underestimate the data points compared to the MDA results of 
Ref. |7| and this effect should be due to the use of single-folding method (see 
further discussion below). 

In the energy bin centered at E x = 10.9 MeV for the 240 MeV data (or 
E x = 10.5 MeV for the 386 MeV data) the E2 strengths deduced from both 
data sets are quite close to each other, exhausting about 19 to 23% of the E2 
EWSR (see Tabled]). In addition, similar E0 and El strengths (from 3 to 
4 % of the corresponding EWSR) were also deduced from both experiments. 
The only difference is a significantly larger E3 strength deduced by the MDA 
of the 240 MeV data for this energy bin. As there is no (a, a') cross section 
measured at E^ = 240 MeV available to us for this energy bin, we compare 
in Fig. [6] the (a, a') data at E^ = 386 MeV measured for the energy bin 
centered at E x = 10.5 MeV |39j with the DFM + DWBA results given by the 
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isoscalar EL strengths taken from Refs. [7|, [8|, [9( . One can see in lower panel 



of Fig. M that the DFM + DWBA resu 
from the MDA of the 240 MeV data 



ts based on the EL strengths taken 
7j are much closer to the data points 



compared with similar results for the (a, a') cross sections at E\ ah = 386 
MeV shown in lower panels of Figs. [3] and 0} This could well indicate a 
much weaker contribution by the high multipoles to the (a, a') cross section 
at the excitation energy around 10 MeV. Like the results obtained above 
for the excitation energies around 14 MeV, the DFM + DWBA calculation 
based on the EL strengths taken from the MDA of the 386 MeV data [s, ^] 
underestimates the data over the whole angular range (see upper panels of 
Figs. El HI and [6]) including the smallest angles where high multipoles are not 
expected to play a major role. 

3.4- Isoscalar EL strengths near the main ISGDR peak 

The ISGDR in 208 Pb has been observed in both the (a, a') experiments at 
240 MeV [?] and 386 MeV jj. The isoscalar El strength distribution has 
been shown clearly by these two experiments to split into two parts: a weak 
low-energy peak centered at E x « 13 MeV and the main, broad high-energy 
peak at E x ~ 22.5 MeV. In contrast to the ISGMR case, the ISGDR peaks 
observed in these two measurements are quite close to each other, except for 
some difference in the width deduced for the low-energy El peak. Using the 



empirical formulas 40j, |41[ relating the high-energy ISGDR peak in 208 Pb 
and nuclear matter incompressibility we obtain « 210 MeV which 
is smaller than that deduced from the ISGMR data by about 20 MeV. 

The inelastic a+ 208 Pb scattering data at E\^ = 240 MeV measured for 
the energy bin centered at E x = 22.5 MeV (3] and DFM + DWBA results 
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given by the CM transition densities scaled to the isoscalar EL strengths 
taken from Refs. are shown in Fig. [7J After substraction of the 

continuum contribution, the MDA [7| implied that the measured (a, a') cross 
section (see Fig. [7J contains mainly the isoscalar El and E3 components 
which exhaust, respectively, about 8% and 6% of the corresponding EWSR 
(see Table |2]). The DFM + DWBA calculation based exactly on these El 
and E3 transition strengths accounts reasonably for the data (lower panel 
of Fig. [7j), with some underestimation of the data points at large angles 
(due perhaps to the contribution from the isovector El mode, see Fig. 5 of 
Ref. [7|). In contrast to the MDA results for the 240 MeV data, in addition 
to similar isoscalar El and E3 strengths around the peak E x = 22.5 MeV, 



the MDA of the 386 MeV data 



has found significant contributions 



from the isoscalar E0 and E2 strengths which exhaust, respectively, about 
2% and 10% of the corresponding EWSR (Table EJ). The DFM + DWBA 
calculation based on the isoscalar EL strengths given by the MDA of the 
386 MeV data also describes reasonably the 240 MeV data for the peak 
E x = 22.5 MeV (upper panel of Fig. [7J. The high-energy tails of the E0 
and El strength distributions given by the MDA of the 386 MeV data give 
rise to an enhancement of the DWBA cross section at the forward angles as 
shown in upper panel of Fig. [7J In Fig. [8] the 386 MeV data measured for 
the 1 MeV bin centered at E x = 22.5 MeV are compared with the DFM 
+ DWBA results given by the CM transition densities fll3l) - (TT4|) based on 
the same isoscalar EL transition strengths as those used in Fig. [7J Similar 
to the DFM + DWBA results shown in Figs. [3] and H the DFM + DWBA 
results for the peak E x = 22.5 MeV strongly underestimate the 386 MeV data 
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over the whole angular range. However, the gap between the calculated and 
measured cross sections becomes significantly larger in this case which shows 
the important contributions by the excitation modes of higher multipoles 
(L > 3) at energies above 20 MeV. 

Given a sizable difference between the DFM + DWBA results and the 
MDA results of Ref. 7j found for the energy bin centered at the ISGMR peak 
shown in Fig. HI it is necessary to check this effect also for the energy bins 
centered at the ISGQR and ISGDR peaks. In Fig. U the DFM + DWBA 
results for the (a, a') cross sections at = 386 MeV (in the 1-MeV bins 
centered at E x = 10.5 and 22.5 MeV) are compared with the corresponding 
MDA results by Uchida et al. [s[ 9] that are based on the same EL strengths. 
Although in logarithmic scale these two sets of calculated DWBA cross sec- 
tions look similar in shape and strength, the MDA cross sections are larger 
than those given by the DFM + DWBA calculation by about 40 ~ 60% over 
the whole angular range. For the energy bin centered at the ISGQR peak, 
this difference is up to ~ 180% at the most forward angles and is due mainly 
to the difference in the E2 cross sections. Since the Coulomb contribution of 
the inelastic scattering FF is quite strong at the forward angles for the E2 
excitation mode, such an unusually large difference in the E2 cross sections 
could be due to the different treatments of the Coulomb inelastic scattering 
FF in the two approaches. Namely, the Coulomb FF is evaluated in our DFM 
approach microscopically by double-folding the proton parts of the 4 He g.s. 
density and 208 Pb transition density with the Coulomb interaction, with both 
the direct and exchange terms calculated by a method similar to that used for 
the nuclear FF, whereas the widely used ansatz for the Coulomb FF in the 
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MDA of (a, a') data is to assume a simple macroscopic model-independent 
formula 1 261 containing the electric transition rate B(EL) of the considered 
state. To probe this effect, we have made the single-folding calculation [15| of 
the nuclear FF for these cases, using exactly the same effective aN and CM 
transition densities as those used in the MDA by Uchida et al. jsl, Q]. The 
single-folded nuclear FF were then used with the same microscopic Coulomb 
FF as that used in the DFM + DWBA calculation to calculate the (a, a') 
cross sections in the energy bins centered at 10.5 and 22.5 MeV and the 
results are plotted in Fig. [9] as dash-dotted curves. One can see that the 
large difference at forward angles in the cross sections given by the single- 
and double-folding methods is reduced significantly, and at large angles the 
cross sections given by the single-folded FF are very close to the MDA cross 
sections. The results shown in Fig. [9] stress again the need to use the accurate 
DFM in the MDA of the (a, a') data to deduce the realistic EL transition 
strengths. 

In conclusion, the DFM + DWBA analysis of the inelastic a+ 208 Pb scat- 
tering data at = 240 MeV measured in the energy bins centered at the 
peaks of the ISGMR, ISGQR and ISGDR in 208 Pb, using the CM nuclear 
transition densities for the EL excitation with L ^ 3, agree qualitatively 
with the original MDA of these data [7|. Given a sizable difference in the 
(a, a') cross sections obtained with the single- and double- folding methods 
for the inelastic scattering FF, the uncertainties in the EL strengths deduced 
from the MDA of the considered data 7|, |8|, |9J might be significantly larger. 
Similar DFM + DWBA analysis of the inelastic a+ 208 Pb scattering data 
at E\ ah = 386 MeV strongly underestimates the data points measured in 
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about the same energy bins and, thus, indicates a strong contribution by the 
high-mult ipole (L > 3) excitation modes. However, the gap between the cal- 
culated and measured cross sections is quite different in the excitation energy 
regions around 10 MeV and above 20 MeV. This result shows that the max- 
imum angular momentum L max taken into account in the MDA series (12171) 
seems to be energy dependent if the background due to the high-multipole 
excitation modes is not explicitly subtracted which is the case for the 386 
MeV data [JSE]- 

4. Results obtained with the microscopic RPA transition densities 

Although there is no consensus whether microscopic models like RPA can 



provide reliab 
in the past 31 



e in put s for the nuclear transition densities, it has been shown 
42j, |43j that for the low-lying excited states of dominating one- 



phonon structure, the RPA transition densities can be successfully used in the 
folding model analysis. The DFM + DWBA description of inelastic a+ 208 Pb 
scattering data measured for the 3]~ state of 208 Pb given by the RPA nuclear 
transition density (see lower panel of Fig. [1]) is again a convincing example. Of 
course, there are cases in which RPA is known to have drawbacks, like in the 
case of low-lying states with a strong anharmonic mixture of the two-phonon 



structure [42| . Isoscalar giant resonances, as already discussed above, should 
be a good test ground for the RPA wave functions because RPA has been 
claimed over the years to be a proper theory to describe those resonances. 
However, these qualitative arguments are often invoked but in fact we are 
not aware of the conclusive evidences showing that microscopic RPA provides 
accurate transition densities for inelastic scattering calculations (cf., in this 
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respect, Refs. [44J). Also, this is probably the first combination of a fully 
self-consistent RPA approach and an advanced microscopic double-folding 
model, and consequently in the present context the question of the accuracy 
of the RPA transition densities can be addressed carefully in more detail. 

It is impossible to reproduce the full experimental width of a resonance 
state within the RPA approach. Although the width caused by fragmenta- 
tion of the resonance strength (the so-called Landau damping) can be ac- 
counted for within the RPA, and the escape width can also be accounted 
for if continuum-RPA is performed, the spreading width (which is by far the 
most relevant in heavy nuclei) cannot come out from RPA. To have a direct 
quantitative comparison of the RPA solution with the observed EL strength 
distribution, we found it necessary to perform some averaging [45( of the 
total RPA strength ([19]) over the excitation energy E x as follows 

{Sr\E x )) = J2 Sr\E v )f{E x - E„), (21) 

V 

n 

where v labels the RPA (isoscalar EL) states and a Lorentzian [45| is used 
as the averaging function f(E — E'). In each case, the averaging width A is 
chosen so that the averaged RPA strength 

pE x +0.5 MeV 

/ (Sr A (E))dE (22) 

JE x -0.5 MeV 

within the 1-MeV energy bin centered at an excitation energy E x is smooth 
enough to be compared with that deduced from the MDA analysis of (a, a') 
data. In the case of ISGMR, a strongly collective resonance RPA state was 
found at E x w 14.2 MeV which is quite close to the experimental ISGMR 
peak. This RPA state is accompanied by several non-collective RPA states 
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on either sides of the peak (see upper panel of Fig. [TU]) and the whole set of 
monopole RPA states exhausts about 99.5% of the EO EWSR. The averaging 
width A si 3 MeV, which is close to the observed ISGMR width [7|, was 
chosen to smooth the strength of the RPA resonance state over the excitation 
energy. The distribution of averaged RPA strength agrees reasonably with 
those deduced from the MDA analyses of the (a, a') data at 240 MeV 7| and 



386 MeV 



39J as shown in the lower panel of Fig. [TUJ 



For the ISGQR, most of the microscopic RPA calculations predict a strong 
E2 resonance at somewhat higher excitation energy (E x « 12.5 MeV in 
our case) compared to the experimental ISGQR peak around 10.3 MeV, 
in keeping with the low effective mass associated with effective Skyrme or 
Gogny interactions. The lowest 2f state is predicted by the RPA at an 
excitation energy of E x « 5.1 MeV which is also higher than the experimental 
value of about 4.09 MeV. To have comparable E2 strengths at the ISGQR 
peak observed in (a, a') experiment, we have shifted all isoscalar quadrupole 
RPA states (which exhaust 99.3% of the E2 EWSR) down by 2 MeV in the 
excitation energy and the resulted RPA spectrum is shown in upper panel of 
Fig. [TTJ By using the same averaging width A « 3 MeV which is also close 
to the observed ISGQR width the distribution of averaged RPA strength 
agrees well with those deduced from the MDA analyses of the ( cy, ol ) data at 
both energies (see lower panel of Fig. [TT|) . While the averaged RPA strength 
and that given by the MDA of 240 MeV data [7] are concentrated mainly 
near the ISGQR peak and slightly spread over E x < 21 MeV, the E2 strength 
given by the MDA of 386 MeV data [s[ is broadly spread up to much 
higher excitation energies. 
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Concerning the ISGDR excitation in Pb, most of the microscopic struc 



ture calculations 
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47, 



48j predict the main (high-energy) peak of the 



ISGDR at E x pa 24.5 ~ 25.5 MeV which is somewhat higher than that 



a la 



39j. Like the 



(E x pa 22.5 MeV) observed in the (a, a') experiments 
earlier RPA results obtained with the SLy4 interaction |46j], the present RPA 
calculation using the SLy5 interaction predicts the high-energy peak of the 
ISGDR at E x pa 25 MeV. In a similar manner, we have shifted all the isoscalar 
dipole RPA states (which exhaust 86.5% of the El EWSR) down by 3 MeV 
in excitation energy and the shifted spectrum of dipole RPA states is shown 
in the upper panel of Fig. [12j To have a better resolution of the averaged 
dipole strength of numerous RPA states found in the resonance region, we 
have used a finer averaging width of A pa 1 MeV in this case and the distri- 
bution of averaged RPA strength agrees reasonably with the observed dipole 
strength distribution at the main ISGDR peak (see lower panel of Fig. [T2|) . 
However, as it can be seen in the lower panel of Fig. [T2J the low-energy IS 
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39) 



GDR strength observed at E x pa 13 MeV in both (a, a') experiments 
is not reproduced by the present RPA calculation. Whether this is an ex- 
perimental or theoretical problem, is an open question. Certainly in most of 
the microscopic RPA calculations the low-energy ISGDR strength is less col- 
lective than the high-energy part (see, e.g., the discussion in Ref. 46] ) and, 
as such, less amenable to a RPA description. Suggestions that this strength 
corresponds to toroidal motion have been put forward 49j and if this were 
the case, the capability of microscopic functionals to describe such exotic 
mode, and the relationship with (a, a') cross sections is even less clear. 



Given the energy distributions of the EO, El and E2 strengths reasonably 
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described by the averaged RPA strengths as shown above, it is natural to 
expect that the DFM + DWBA calculation using a proper input of the RPA 
transition densities should also deliver a good description of the measured 
(a, a') cross sections. However, to compare the DFM + DWBA results with 
the (a, a') cross section measured for a given energy bin, one needs to combine 
properly the transition densities of all RPA states in this energy bin into 
a total RPA transition density which can be used as input of the DFM 
calculation. In general, if the number of RPA solutions for a given EL 
transition in the energy bin dE centered at E x is N, then the total RPA 
transition density associated with the bin dE should be defined as 

N N N 

$pl PA (r) = M Lu 5p Lu (r); with S* PA = S^ A = \M Lu \ 2 . (23) 

Here Ml„ is the RPA transition moment (fl9l) of the RPA state \Lv > and 
8pLv( r ) is the corresponding RPA transition density (ITS]) . It is natural to 
choose the transition moment Mi u as the averaging weight for the RPA tran- 
sition density 5pi v (r), so that the RPA transition density fl23|) preserves the 
total transition strength S^f in the considered energy bin as predicted by 
the RPA calculation. The only question now is whether the total RPA tran- 
sition density should be scaled to reproduce the averaged RPA strength (1221) 
in this energy bin or should it be kept unchanged as given by the weighted 
sum f !23p of all RPA transition densities. We have found, however, that the 
first procedure is reasonable only for the transition densities 5pL„(r) of the 
strongest RPA states with L = and 2 in the energy bins around the ISGMR 
and ISGQR peaks, respectively. For the RPA states of other multipolarities 
(L — 1, 3) in the same energy bins the RPA spectrum consists only of a few 
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non-collective states and such a scaling procedure can lead to very unrealis- 
tic shapes of SpLu{r) with L = 1 and 3 which strongly distort the calculated 
(a, a') cross section for the ISGMR and ISGQR peaks. We have used, there- 
fore, the total RPA transition density as given by the weighted sum f l23l) 
of all RPA transition densities in the considered energy bin for the DFM + 
DWBA calculation of the (a, a') cross section. 

To compare our DFM + DWBA results with the inelastic a+ 208 Pb scat- 
tering data at Ei ab = 240 MeV [7J, we have generated the total RPA transi- 
tion densities (123]) for the three 640-keV energy bins centered at E x = 10.3, 
14.1 and 22.5 MeV, respectively. The total RPA transition density and the 
two strongest individual RPA transition densities (fl~8j) in each energy bin 
are compared in Fig. [T3] with the total collective model transition density 
(CEID-CLSD based on EL strengths given by the MDA of the 240 MeV data [7|. 
While the radial shape of the total RPA transition densities in the energy 
bins centered at the ISGQR and ISGMR peaks agrees more or less with that 
of the CM densities, the EL strengths given by the RPA transition densities 
are much stronger than those given by the CM densities. Such an effect is not 
unexpected because the predicted E0 and E2 strengths near the ISGQR and 
ISQMR peaks are concentrated in just a few discrete RPA states (see upper 
panels of Figs. [TU1 and [TTj) . In contrary, the RPA strength for the isoscalar 
dipole excitation are distributed over many weakly excited El states, and 
the total RPA dipole transition density given by the weighted sum (1231) of 
15 RPA states in the 640-keV bin centered at E x = 22.5 MeV has a slightly 
weaker El strength compared to that of the CM transition density. In other 
words, the RPA calculations of the E0 and E2 strength do not predict frag- 
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mentation, at variance with the El case where the fragmentation caused by 
Landau damping is quite large. The EL transition strengths given by the 
total RPA transition densities fl23l) in the three 640-keV energy bins centered 
at the ISGMR, ISGQR and ISGDR peaks are given in Table [31 

The DFM + DWBA description of the (a, a') data at 240 MeV given 
by the RPA transition densities are shown in Fig. [TH Although these DFM 
+ DWBA results agree well with the data, in about the same way as the 
results given by the CM transition densities, the EL strengths associated 
with the RPA transition densities in each energy bin are quite different from 
those given by the CM transition densities (see Table [2]). In energy bins 
centered at the ISGMR and ISGDR peaks the EL transition strengths with 
L = and 2, respectively, are much more dominant compared to those 
deduced from the MDA of the (a, a') data. For example, only a single RPA 
transition density (fT8l) of the strongest (discrete) 2 + state at the ISGQR 
peak accounts perfectly for the data measured for the energy bin centered at 
E x = 10.3 MeV. However, this result should not give the wrong impression 
that the EL strengths from other multipoles are not as significant as given 
by the MDA of 240 MeV data [7]. In reality, the isoscalar giant resonances 
under study are not discrete states as predicted by the RPA but widely 
fragmented over the excitation energy, having widths r pa 3 ~ 4 MeV. The 
physical origin of the observed EL strength fragmentation of the isoscalar 
giant resonances, in particular their spreading widths, can be described by 
invoking the anharmonic effects beyond RPA caused by, e.g., the coupling 



with 2p-2h type configurations (see Refs. 50|, |5lJ for extensive reviews). 



The inclusion of the anharmonic effects is expected not only to redistribute 
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the strength of the resonance RPA state over the energy like the averaging 
procedure fl2Tl) . but also to pull down the predicted excitation energy to 
a lower value. We recall again here that the RPA spectrum of 2 + states 
has been shifted down by 2 MeV in energy to have the strongest 2 + state 
near the observed ISGQR peak (see upper panel of Fig. [LTJ); this shift is 
precisely associated with 2p-2h coupling or in other words, as stated above, 
to the renormalization of the effective mass m*. Nevertheless, such a perfect 
agreement of the DWBA cross section predicted by the total RPA transition 
density with the data measured for the ISGQR peak (see upper panel of 
Fig. [T4"j) without any fine-tuning of the FF strength is very encouraging. 
This result confirms the realistic shape of the RPA transition density shown 
in Fig. [T3] and stresses once more the strong predicting power of the RPA 
approach. 

For the ISGMR, the DFM + DWBA calculation using a single transition 
density (ITS]) for the collective RPA state near the ISGMR peak predicts 
values of the (a, a') cross section quite close to the data in magnitude but 
with a deep oscillation pattern that can be smoothed out only by adding 
the contributions from other multipoles. Such a behavior of the (a, a') cross 
section calculated for the ISGMR excitation in 208 Pb has been seen in the 



earlier folding model studies 
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52| . Although the non-zero multipole EL 
transition strengths predicted by the RPA for the energy bin centered at 
E x = 14.1 MeV are somewhat weaker than those given by the MDA of 
240 MeV data (see Tables [2] and [3]) , their contributions are still essential in 
smoothing out the deep oscillation of the DWBA cross section predicted by 
the E0 transition. Given the quite accurate RPA description of both the 
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EO strength distribution (upper panel of Fig. [TU]) and of the (a, a') cross 
section measured for the ISGMR peak (middle panel of Fig. [T4"j) without any 
readjustment of the model parameters and shift of the excitation energy, we 
conclude that the RPA is indeed a reliable theoretical approach to study the 
GMR excitations in nuclei. 

For the ISGDR, the DFM + DWBA calculation using the total RPA 
transition densities (fT8|) for L < 3 in the energy bin centered at E x = 22.5 
MeV reproduces very well the (a, a') angular distribution measured at 240 
MeV (see lower panel of Fig. [T41) . The main difference between these results 
and the DFM + DWBA results given by the CM transition densities shown 
in Fig. [7| is that the RPA predicts quite a strong E3 strength in this energy 
bin (~ 10.8%) compared to that deduced from the MDA of 240 MeV data 
(~ 5.9%), while the El strength predicted by the RPA is about 60% weaker 
than the MDA value. As a result, the E3 transition turned out to give a 
dominant contribution to the DWBA cross section for the energy bin centered 
at E x = 22.5 MeV. 

In conclusion, our study of the ISGMR, ISGQR and ISGDR strength 
distributions and inelastic scattering (a, a') angular distributions measured 
at the corresponding resonance peaks using the RPA transition densities has 
shown that the RPA approach can be successfully used to describe not only 
the EL strength distribution and integral properties of the giant resonances, 
but also the (double) differential cPa/dQdE cross sections. Although there 
are some differences in the EL strengths of non- collective states, the ISGMR, 
ISGQR and ISGDR strength distributions predicted by the averaged RPA 
results agree reasonably good with those deduced from the MDA of the (a, a') 
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data (see Tableland Figs. [TUJ [TT1 and IT5T) . The EL transition strengths as 
well as the location of the ISGQR and ISGDR peaks originally predicted 
by the RPA are significantly higher than those given by the MDA of the 
considered (a, a 1 ) data. This effect can be qualitatively explained by the 
lack of the (beyond RPA) anharmonic contributions from, e.g., the 2p-2h 
coupling. 

Therefore, it would be of further interest to perform the same DFM + 
DWBA analysis including effects beyond the RPA since a self-consistent mi- 
croscopic calculation beyond RPA coupled to an accurate reaction framework 
does not exist yet, to our knowledge. The effects, which have been introduced 
ad hoc in the present work via the shift of the RPA mean peak(s) downwards 
and spreading (|2ip - (|22|) of the discrete RPA strength, are exactly those in- 
cluded microscopically into a second RPA, or RPA plus phonon coupling 



(RPA-PC) calculations (see also Ref. 



). In this sense, the present results 



are encouraging. In particular, the effects of 2p-2h coupling on the ISGMR 
are expected to be small due to well-known cancellation effects and this goes 
along with the fact that our present simple RPA gives a very good DFM + 
DWBA description of the monopole excitation. 

5. Summary 

The generalized double-folding approach of Refs. Q, [r| has been further 
developed for the microscopic study of isoscalar giant resonances induced by 
inelastic a-nucleus scattering, using the nuclear transition densities given by 
both the collective model and microscopic self-consistent RPA calculation. 

n 

Although the single-folding approach [15] has been often used to compute 
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the a-nucleus potential in the MDA of the (a, a') data, we have shown in the 
present work a significant difference in the (a, a') cross sections given by the 
single- and double-folding methods. The DFM has also been suggested earlier 



m 



26| as a more accurate approach to obtain realistic results for the inelastic 



a-nucleus and HI scattering, in particular, the EWSR fractions exhausted 
by different EL excitation modes. Therefore, the present combination of the 
DFM and DWBA approaches should be a reliable alternative method to be 
used in the MDA of future (a, a 1 ) measurements. 

A reasonable DFM + DWBA description of the 240 MeV (a, a') data in 
energy bins centered at the ISGMR, ISGQR and ISGDR peaks in 208 Pb has 
been obtained with the CM nuclear transition densities built upon the EL 
strengths given by the MDA of these data. Our similar study of the inelastic 
a+ 208 Pb scattering at E\ a \ } = 386 MeV strongly underestimates the data 
points measured in the same energy bins and, thus, shows quite a strong 
contribution by the excitation modes of higher multipole (L > 3). The 
contribution by higher multipole excitation modes was found significantly 
different in the two energy regions around 10 MeV and above 20 MeV, which 
suggests that the maximum angular momentum L max taken into account in 
the MDA series ( 12"U|) should be energy dependent. 

The present DFM + DWBA method has also provided an accurate di- 
rect link between the discrete RPA approach (which was used in the past to 
describe mainly the integral properties of the giant resonances) and the exper- 
imental double differential (a, a') cross sections measured for the resonance 
peaks. Given high-precision (a, a') data for the isoscalar giant resonances 
under study, our method can be used in the future to probe the accuracy of 
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the microscopic prediction by different structure models for the energy dis- 
tribution of the Sl(E) strength as well as the inelastic scattering d 2 a/dVLdE 
cross section. In such a connection, we would like to emphasize again that 
the latter is very sensitive to the interference of different EL contributions. 

Acknowledgments 

We are indebted to Prof. D.H. Youngblood, Prof. T. Kawabata, Dr. X. 
Chen and Dr. M. Uchida for their helpful communications on the (a, a') 
data under study. The present research has been supported by National 
Foundation for Scientific & Technological Development (NAFOSTED) un- 
der Project No 103.04.07.09. D.C.C. gratefully acknowledges the financial 
support from the Asia Link Programme CN/Asia-Link 008 (94791) during 
his two research stays at the University of Milan where this work has been 
initiated. Last but not least, one of us (D.T.K.) is grateful to Ray Satchler 
for his numerous discussions during the years of collaboration, which have 
laid the foundation for this research study. 

References 

[1] M.N. Harakeh, A. van der Woude, Giant Resonances: Fundamental 
High- Frequency Modes of Nuclear Excitation, Clarendon Press, Oxford, 
2001; P.F. Bortignon, A. Bracco, R.A. Broglia, Giant Resonances; Nu- 
clear Structure at Finite Temperature, Harwood Academic, New York, 
1998. 

[2] S. Shlomo, V.M. Kolomietz, G. Colo, Eur. Phys. J. A 30 (2006) 23-30. 



40 



[3] G. Colo, Physics of Particles and Nuclei 39 (2008) 286-326. 

[4] M.N. Harakeh, K. van der Borg, T. Ishimatsu, H.P. Morsch, A. van der 
Woude, and F.E. Bertrand, Phys. Rev. Lett. 38 (1977) 676-679. 

[5] M.N. Harakeh, B. van Heyst, K. van der Borg, A. van der Woude, Nucl. 
Phys. A 327 (1979) 373-396. 

[6] D.H. Youngblood, CM. Rozsa, J.M. Moss, D.R. Brown, J.D. Bronson, 
Phys. Rev. Lett. 39 (1977) 1188-1191. 

[7] D.H. Youngblood, Y.W. Lui, H.L. Clark, B. John, Y. Tokimoto, X. 
Chen, Phys. Rev. C 69 (2004) 034315; D.H. Youngblood, X. Chen, pri- 
vate communication (unpublished). 

[8] M. Uchida, H. Sakaguchi, M. Itoh, M. Yosoi, T. Kawabata, H. Takeda, 
Y. Yasuda, T. Murakami, T. Ishikawa, T. Taki, N. Tsukahara, S. 
Terashima, U. Garg, M. Hedden, B. Kharraja, M. Koss, B. K. Nayak, S. 
Zhu, M. Fujiwara, H. Fujimura, K. Hara, E. Obayashi, H.P. Yoshida, H. 
Akimune, M.N. Harakeh, M. Volkerts, Phys. Lett. B 557 (2003) 12-19. 

[9] M. Uchida, H. Sakaguchi, M. Itoh, M. Yosoi, T. Kawabata, Y. Ya- 
suda, H. Takeda, T. Murakami, S. Terashima, S. Kishi, U. Garg, P. 
Boutachkov, M. Hedden, B. Kharraja, M. Koss, B.K. Nayak, S. Zhu, 
M. Fujiwara, H. Fujimura, H.P. Yoshida, K. Hara, H. Akimune, M.N. 
Harakeh, Phys. Rev. C 69 (2004) 05 1301 (R); T. Kawabata, M. Uchida, 
private communication (unpublished). 

[10] M.N. Harakeh, A.E.L. Dieperink, Phys. Rev. C 23 (1981) 2329-2334. 



41 



[11] N.V. Giai, H. Sagawa, Nucl. Phys. A 371 (1981) 1-18. 

[12] B.F. Davis, U. Garg, W. Reviol, M.N. Harakeh, A. Bacher, G.P. A. 
Berg, C.C. Foster, E.J. Stephenson, Y. Wang, J. Janecke, K. Pham, 
D. Roberts, H. Akimune, M. Fujiwara, J. Lisantti, Phys. Rev. Lett. 79 
(1997) 609-612. 

[13] G.R. Satchler, Direct Nuclear Reactions, Clarendon Press, Oxford, 1983. 

[14] D.T. Khoa, G.R. Satchler, Nucl. Phys. A 668 (2000) 3-41. 

[15] G.R. Satchler, D.T. Khoa, Phys. Rev. C 55 (1997) 285-297. 

[16] D.T. Khoa, G.R. Satchler, W. von Oertzen, Phys. Rev. C 56 (1997) 
954-969. 

[17] D.T. Khoa, W. von Oertzen, H.G. Bohlen, S. Ohkubo, J. Phys. G 34 
(2007) R111-R164. 

[18] J.P Jeukenne, A. Lejeune, C. Mahaux, Phys. Rev. C 16 (1977) 80-96. 

[19] D.T. Khoa, D.C. Cuong, Phys. Lett. B 660 (2008) 331-338. 

[20] N. Anantaraman, H. Toki, G.F. Bertsch, Nucl. Phys. A 398 (1983) 269- 
278. 

[21] A. Ingemarsson, J. Nyberg, P.U. Renberg, O. Sundberg, R.F. Carlson, 
A.J. Cox, A. Auce, R. Johansson, G. Tibell, D.T. Khoa, R. E. Warner, 
Nucl. Phys. A 676 (2000) 3-31. 



42 



[22] B. Bonin, N. Alamanos, B. Berthier, G. Bruge, H. Faraggi, J.C. Lu- 
gol, W. Mittig, L. Papineau, A.I. Yavin, J. Arvieux, L. Farvacque, M. 
Buenerd, W. Bauhoff, Nucl. Phys. A 445 (1985) 381-407. 

[23] D.T. Khoa, Phys. Rev. C 63 (2001) 034007. 

[24] W.G. Love, Nucl. Phys. A 312 (1978) 160-176. 

[25] F. Carstoiu, M. Lassaut, Nucl. Phys. A 597 (1996) 269-297. 

[26] D.J. Horen, J.R. Beene, G.R. Satchler, Phys. Rev. C 52 (1995) 1554- 
1564. 

[27] A. Bohr, B.R. Mottelson, Nuclear Structure, W.A. Benjamin, Reading, 
1975, vol. II. 

[28] M.E. Farid, G.R. Satchler, Nucl. Phys. A 438 (1985) 525-535. 

[29] H. Uberall, Electron Scattering from Complex Nuclei, Academic Press, 
New York, 1971, vol. B. 

[30] A. van der Woude, Nucl. Phys. A 649 (1999) 97c-103c. 

[31] D.T. Khoa, E. Khan, G. Colo, N.V. Giai, Nucl. Phys. A 706 (2002) 
61-84. 

[32] J. Guillot, S. Gales, D. Beaumel, S. Fortier, E. Rich, N. Van Giai, G. 
Colo, A.M. van der Berg, S. Brandenburg, B. Davids, M.N. Harakeh, M. 
Hunyadi, M. de Huu, S.Y. van der Werf, H.J. Wortche, C. Baumer, D. 
Frekers, E.W. Grewe, P. Haefner, B.C. Junk, M. Fujiwara, Phys. Rev. 
C 73 (2006) 014616. 

43 



[33] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, R. Schaeffer, Nucl. Phys. 
A 635 (1998) 231-256. 

[34] G. Colo, PF. Bortignon, S. Fracasso, N. V. Giai, Nucl. Phys. A 788 
(2007) 137c-181c. 

[35] G.R. Satchler, W.G. Love, Phys. Rep. 55 (1979) 183-254. 

[36] T. Kibedi, R.H. Spear, At. Data and Nucl. Data Tables 80 (2002) 35-82. 

[37] G. Colo, N. Van Giai, J. Meyer, K. Bennaceur, P. Bonche, Phys. Rev. 
C 70 (2004) 024307. 

[38] D. Vretenar, T. Niksic, P. Ring, Phys. Rev. C 68 (2003) 024310. 

[39] M. Uchida, PhD Thesis (University of Kyoto, 2003). 

[40] S. Stringari, Phys. Lett. 108 B (1982) 232-236. 

[41] J.P. Blaizot, J.F. Berger, J. Decharge, M. Girod, Nucl. Phys. A 591 

(1995) 435-457. 

[42] D.T. Khoa, I.N. Kuchtina, V.Yu. Ponomarev, Sov. J. Nucl. Phys. 44 
(1986) 585-590. 

[43] D.J. Horen, G.R. Satchler, S.A. Fayans, E.L. Trykov, Nucl. Phys. A 600 

(1996) 193-235. 

[44] A. Kolomiets, O. Pochivalov, S. Shlomo, Phys. Rev. C 61 (2000) 034312. 

[45] H. Feshbach, Theoretical Nuclear Physics - Nuclear Reactions, Wiley, 
New York, 1992. 

44 



[46] G. Colo, N.V. Giai, P.F. Bortignon, M.R. Quaglia, Phys. Lett. B 485 
(2000) 362-366. 

[47] D. Vretenar, A. Wandelt, P. Ring, Phys. Lett. B 487 (2000) 334-340. 
[48] J. Piekarewicz, Phys. Rev. C 64 (2001) 024307. 

[49] S.I. Bastrukov, S. Misicu, V.I. Sushkov, Nucl. Phys. A 562 (1993) 191; S. 
Misicu, Phys. Rev. C 73 (2006) 024301; D. Vretenar, N. Paar, P. Ring, 
T. Niksic, Phys. Rev. C 65 (2002) 021301(R). 

[50] G.F. Bertsch, P.F. Bortignon, R.A. Broglia, Rev. Mod. Phys. 55 (1983) 
287-314. 

[51] S. Drozdz, S. Nishizaki, J. Speth, J. Wambach, Phys. Rep. 197 (1990) 
1-65. 

[52] F.E. Bertrand, G.R. Satchler, D.J. Horen, J.R. Wu, A.D. Bacher, G.T. 
Emery, W.P. Jones, D.W. Miller, A. van der Woude, Phys. Rev. C 22 
(1980) 1832-1847. 

[53] V. Tselyaev, J. Speth, S. Krewald, E. Litvinova, S. Kamerdzhiev, 
N. Lyutorovich, A. Adveenkov, F. Grummer, Phys. Rev. C79 (2009) 
034309. 



45 



Figure captions 



Fig.l: Upper panel: OM description of the elastic a+ Pb scattering 



data at 240 MeV 



7j and 386 MeV 8|, |9] given by the unrenormalized (dash 



7 



ower panel: 
, b|, b| for 3^ 



curves) and renormalized complex folded OP (solid curves). 
DWBA descriptions of the inelastic a+ 208 Pb scattering data 
state of 208 Pb (E x = 2.61 MeV) given by the collective model (dash curves) 
and RPA (solid curves) nuclear transition densities. 

Fig. 2: Inelastic a+ 208 Pb scattering data at E^b = 240 MeV measured for 
the energy bin centered at E x = 14.1 MeV jjl] in comparison with the DFM 
+ DWBA results given by the CM transition densities based on the isoscalar 
EL strengths taken from Refs. [?] (lower panel) and [8f (upper panel). See 
details in the text and Table [2J 

Fig. 3: Inelastic ct+ 208 Pb scattering data at E ]?h = 386 MeV measured 
for the energy bin centered at E x = 13.5 MeV 39( in comparison with the 
DFM + DWBA results given by the CM transition densities based on the 



7J (lower panel) and 



isoscalar EL strengths taken from Refs. 
panel). See details in the text and Table [21 

Fig. 4: Inelastic a+ 208 Pb scattering data at E^h 



aa 



(upper 



386 MeV measured 



for the energy bin centered at E x = 14.5 MeV 



391 ] in comparison with the 
DWBA results given by the isoscalar EL strengths taken from Refs. 8|, |9|. 
Upper panel: the present DFM + DWBA calculation; lower panel: the MDA 
results by Uchida et al. js, 9]. 

Fig. 5: The same as Fig. [2] but for the energy bin centered at E x = 10.3 
MeV. 
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Fig. 6: The same as Fig. [3] but for the energy bin centered at E x = 10.5 
MeV. 

Fig. 7: The same as Fig. [2] but for the energy bin centered at E x = 22.5 
MeV. 

Fig. 8: The same as Fig. [3] but for the energy bin centered at E x = 22.5 
MeV. 

Fig. 9: Inelastic a+ 208 Pb scattering data at = 386 MeV measured 
for the energy bins centered at E x = 10.5 and 22.5 MeV [39J] in comparison 
with the DWBA results given by the isoscalar EL strengths taken from 
Refs. 0,3- Solid curves: the present DFM + DWBA calculation; dashed 
curves: the MDA results by Uchida et al. 8|, [9j; dashed-dotted curves: the 
same as dashed curves but using microscopic Coulomb FF from the present 
DFM calculation. 

Fig. 10: Isoscalar E0 strength distributions deduced from the MDA anal- 
yses of inelastic a+ 208 Pb scattering data at 240 MeV by Youngblood et al. 
7J and 386 MeV by Uchida et al. |8j, |39j in comparison with the RPA results. 
See details in text. 

Fig. 11: The same as Fig. [10] but for the isoscalar E2 strength distribu- 
tions. 

Fig. 12: The same as Fig. [10] but for the isoscalar El strength distribu- 
tions. 

Fig. 13: Total RPA transition density (1231) and transition densities of the 
two strongest RPA states in the 640-keV energy bins centered at E x = 10.3, 
14.1 and 22.5 MeV, respectively. The corresponding collective model tran- 
sition densities were built upon the EL strengths given by the MDA of 240 
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a. 



MeV data [7|. The quoted percentages are the exhausted fractions of the 
isoscalar EL EWSR. 

Fig. 14: Inelastic a+ 208 Pb scattering data at = 240 MeV measured 
for the 640 keV energy bins centered at E x = 10.3, 14.1 and 22.5 MeV 
respectively 3], in comparison with the DFM + DWBA results obtained 
with the total RPA transition densities (123]) which give the fractions of the 
isoscalar EL EWSR shown in Table [31 
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Table 1: Parameters of the complex density dependence of the CDM3Y6 interaction (j3])-(j4| 
used to calculate the OP and inelastic scattering FF for the elastic and inelastic a+ 208 Pb 
scattering at E\ a b — 240 and 386 MeV. iVj^j) are the renormalization coefficients of the 
real and imaginary OP given by the optical model analysis of elastic scattering data; ctr 
is the calculated total reaction cross section. 
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a ) Experimental total reaction cross section measured at Ej a b = 192 MeV 
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Table 2: Fractions of the isoscalar EL EWSR exhausted in energy bin centered at the 
excitation energy E x deduced from the MDA analyses of inelastic a+ 208 Pb scattering 

q nn 

data at E iah = 240 [7[ and 386 MeV [8J, |g|, and from the averaged RPA results. The 
averaging of RPA transition densities (<RPA>) is discussed below in Sec. SJ 8l are the 
deformation lengths for the CM nuclear transition densities (fl"Tj) - (fT6| based on the EL 
strengths taken from Ref. jjj. 
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16.2 


0.0243 






1" 
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1.85 


0.0016 
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1.48 
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2.22 ± 1.26 


1.29 
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\- 


8.23 ±0.15 


8.67 ±0.59 


7.06 


0.0016 






2+ 


0.0 


10.4 ± 1.70 


0.60 


0.0 






3" 


5.90 ±0.30 


5.42 ±0.95 


4.74 


0.1584 
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Table 3: Fractions of the isoscalar EL EWSR exhausted in the 640- keV energy bin centered 
at the excitation energy E x determined from the RPA transition strengths. 



% EWSR/MeV 


17 


E x = 10.3 MeV 


E x = 14.1 MeV 


E x = 22.5 MeV 


0+ 


1.25 


55.1 


0.06 


1" 


0.32 


2.25 


5.45 


2+ 


61.7 


0.34 


0.06 


3" 


0.69 


0.83 


10.8 
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Figure 1: Upper panel: OM description of the elastic a+ 208 Pb scattering data at 240 
MeV p| and 386 MeV [8j, [9( given by the unrenormalized (dash curves) and renormalized 
complex folded OP (solid curves). Lower panel: DWBA descriptions of the inelastic 
a+ 208p b sca ttering data Q, S, for 3j" state of 208 Pb (E x = 2.61 MeV) given by the 
collective model (dash curves) and RPA (solid curves) nuclear transition densities. 
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Figure 2: Inelastic a+ 208 Pb scattering data at E\ a b = 240 MeV measured for the energy 
bin centered at E x — 14.1 MeV [7[ in comparison with the DFM + DWBA results given 
by the CM transition densities based on the isoscalar EL strengths taken from Refs. Q 
(lower panel) and [8( (upper panel). See details in the text and Table [2] 
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Figure 3: Inelastic a+ 208 Pb scattering data at E\ a b — 386 MeV measured for the energy 
bin centered at E x = 13.5 MeV [39| in comparison with the DFM + DWBA results given 
by the CM transition densities based on the isoscalar EL strengths taken from Refs. Q 
(lower panel) and 8, 9] (upper panel). See details in the text and Table [2] 
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Figure 4: Inelastic a+ 208 Pb scattering data at E\ a b — 386 MeV measured for the energy 
bin centered at E x = 14.5 MeV [39( in comparison with the DWBA results given by the 
isoscalar EL strengths taken from Refs. 0,0]. Upper panel: the present DFM + DWBA 
calculation; lower panel: the MDA results by Uchida et al. 0,0. 
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Figure 5: The same as Fig. [2] but for the energy bin centered at E x = 10.3 MeV. 
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Figure 6: The same as Fig. |3] but for the energy bin centered at E x = 10.5 MeV. 
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Figure 7: The same as Fig. [2] but for the energy bin centered at E x = 22.5 MeV. 



58 



10 
10 2 
10 1 
10° t- 



o 



oo 



o 



> 



LU 

a 

"D 

<M 

■o 



\/ vA.^v^ 



10- r 



10* 
10 1 



10" r 



10- r 



10 



i 1 1 i i i i i i i i i i i 1 1 1 i i 1 1 1 i 



i i 1 1 i in 



208r 



Pb(a,a') 208 Pb*@386 MeV, E =22.5 MeV 



co o 



Trans, strengths by Uchida et al. 



OCDO 



CO 



o 



Co 



o 



Q 



D 



O 



CD 



o 



/ \ 



• \ / \ '•' \ / 
: \ i ,'"\ V \ .---1 

i ii \ ' » . v y 

10 J 1 1 1 1 i V | 1 1 1 1 1 1 \ \ \ 1 1 1 1 1 1 1 1 ft 1 1 1 1 1 1 f i i ' i 1 1 1 1 1 M 1 1 i' | i l i' i 1 1 1 1 1 1 1 1 1 1 1 ; 



Trans, strengths by Youngblood et al. 

co o 



OCDO o o 



L=0 
L=1 
1=2. 
L=3 
Total 




1 1 ■ ■ i ' 



' 







6 8 10 
©c. m .( de 9) 



12 14 



16 



Figure 8: The same as Fig. |3] but for the energy bin centered at E x = 22.5 MeV. 
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Figure 9: Inelastic a+ 208 Pb scattering data at £i a b = 386 MeV measured for the energy 
bins centered at E x = 10.5 and 22.5 MeV 3jJ in comparison with the DWBA results given 
by the isoscalar EL strengths taken from Refs. [g, {J. Solid curves: thepresent DFM + 
DWBA calculation; dashed curves: the MDA results by Uchida et al. [8j, [9(; dashed-dotted 
curves: the same as dashed curves but using microscopic Coulomb FF from the present 
DFM calculation. 
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Figure 11: The same as Fig.[TU]but for the isoscalar E2 strength distributions. 
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Figure 12: The same as Fig. [TU] but for the isoscalar El strength distributions. 
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Figure 13: Total RPA transition density (|23|) and transition densities of the two strongest 
RPA states in the 640-keV energy bins centered at E x = 10.3, 14.1 and 22.5 MeV, re- 
spectively. The corresponding collective model transition densities were built upon the 
EL strengths given by the MDA of 240 MeV data [3]. The quoted percentages are the 
exhausted fractions of the isoscalar EL EWSR. 
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Figure 14: Inelastic a+ 208 Pb scattering data at E\ a ^ = 240 MeV measured for the 640 
keV energy bins centered at E x = 10.3, 14.1 and 22.5 MeV respectively [7[, in comparison 
with the DFM + DWBA results obtained with the total RPA transition densities ([23]) 
which give the fractions of the isoscalar EL EWSR shown in Table [21 
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